Field-induced Coulomb Blockade and Quantum Confinement in Graphene 
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We report an approach to confine the carriers in single-layer graphene, which leads to quantum devices with 
field-induced quantum confinement. We demonstrated that the novel Coulomb-blockade effect evolves under a 
uniform magnetic field perpendicular to the graphene device. Our experimental results show that field-induced 
quantum dots are realized in graphene, and a quantum confinement-deconfinement transition is switched by the 
magnetic field. 



PACS numbers: 72.80.Vp, 73.23.Hk, 73.63.-b, 75.70.Ak 

Graphene consists of a single layer of carbon atoms and 
has provided an attractive stage for studying novel two- 
dimensional electron gases (2DEGs) that exhibit a Dirac- 
particle-type gapless ('massless') linear ('relativistic') energy 
dispersion (TJ-O. Triggered by recent progress in graphene 
physics and technology, immense effort has been devoted to 
fabricate graphene nanodevices such as quantum wires or 
quantum dots (QDs). However, confining massless Dirac 
fermions in graphene is difficult due to Klein tunneling and the 
zero-band-gap electronic structure 13). Therefore, although 
attempts have been made to design graphene devices using 
quantum confinement, they often suffer from severe design 
limitations. They basically consist of small quantum-dot is- 
lands, which confine electrons geometrically, to which nar- 
row graphene-constrictions are connected lISMTTTl. In this case, 
the device performance has been limited due to detailed con- 
striction and edge orientation. Recently, several bottom-up 
processes to fabricate nanoribbons have been developed that 
control the edge structure at the atomic level lfT2lfl4l . But a 
lot of research is still needed to develop fabrication processes 
capable of controlling the atomic- scale edge orientation for a 
given geometry. 

To address this situation, we propose herein a new de- 
vice structure, in which graphene mesoscopic islands are per- 
fectly isolated and metallic contacts are directly deposited 
onto them without constrictions. Such a configuration is free 
from disturbances due to the structural fluctuation of con- 
strictions and allows direct contact to the mesoscopic 2DEG 
system. We present an experimental demonstration of a 
'field-induced' Coulomb-blockade effect and quantum con- 
finement in the graphene device. This confinement is in- 
duced by both a uniform magnetic field perpendicular to the 
graphene sheet and an electrostatic surface-potential formed 
by the metal/graphene junction. Our results indicate that a 
Dirac-type particle is confined, and a controlled quantum- 
confinement-deconfinement transition is realized. 

The graphene samples were prepared by micromechanical 
cleavage of Kish graphite deposited on the surface of a sil- 
icon substrate with a 90 nm oxidized silicon. From these 
samples, we selected a number of graphene flakes by opti- 



cal microscope contrast and Raman spectroscopy measure- 
ments [031 US). Using these techniques, we confirmed that 
the samples used for the present study consisted of a single- 
layer graphene. The isolated-graphene-mesoscopic struc- 
tures were patterned by electron-beam (EB) lithography us- 
ing polymethylmethacrylate (PMMA) resist as an etch mask, 
with O2 reactive ion etching being used to etch away the un- 
protected graphene. Next, EB lithography was again used to 
deposit metal electrodes on the graphene mesoscopic struc- 
tures to fabricate the source-drain contacts. Finally, the sam- 
ples were annealed in Ar and H2 (3 %) at 300 °C for 5 min. 
Moreover, in order to use as a reference of our graphene flake, 
we also made a Hall-bar device on it. From the four-terminal 
measurement, we observed the integer quantum Hall effect of 
a single-layer graphene at low temperatures, and estimated the 
flake mobility to be « 2,500 cm 2 /Vs, the mean free path « 30 
nm at 10 12 cm" 2 hole- and electron-carrier densities. 

Figure 1(a) shows a scanning electron microscopic image 
of a single-layer graphene device. The disk diameter is 550 
nm, and the distance between contacts is L = 200 nm. A 
side-view schematic of the device is shown in Fig. 1(b). A 
highly-/?-doped silicon substrate was used as the back gate, 
with a 90 nm oxidized silicon layer serving as the gate di- 
electric. To measure their electrical properties, the graphene 
devices were mounted in an 3 He cryostat equipped with a su- 
perconducting magnet. A magnetic field (B) was applied per- 
pendicular to the sample chip. All electrical-transport mea- 
surements were carried out at base temperatures of 0.23-1.7 
K. The two-terminal conductance G through the graphene de- 
vice was measured as a function of the gate voltage V g by 
applying a fixed DC bias voltage V S( j. In addition, we mea- 
sured the source-drain current / as a function of both V sd and 
V g and obtained the differential conductance dI/dV S d by nu- 
merically differentiating the 7-V S d curve. Figure 1(c) shows 
the conductance through the graphene mesoscopic device as 
a function of V g at V sd = 100 fiV, B = T and T = 1.7 
K. The conductance curve is asymmetric and V-shaped with 
fluctuations, though we consider that the very low conduc- 
tance is due to poor contact between the metal (titanium) and 
the graphene ifTTl . It has a broad minimum around the Dirac 
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FIG. 1. (color). Configuration and transport characteristics of 
graphene device, (a) Scanning microscope image of the measured 
device with electrode assignment. The scale bar is 100 nm. (b) Side- 
view schematic of the device with an electrode. The axis directions 
are shown. The midpoint between the metal contacts is defined as 
x = 0, and the center of the width of the metal segment is defined as 
y = 0. (c) Conductance (measured at V sd = 100 /uV) as a function of 
V g at T = 1.7 K and 5 = 01 (d) Conductance as a function of V g 
under magnetic fields applied perpendicular to the graphene device 
at T = 1.7 K and V sd = 100 //V. For clarity, each trace is shifted in 
100 nS steps from the bottom to the top. 



point at V g = VDirac ~ +4 V. This asymmetry and shift in the 
gate voltage of the Dirac-point energy have been reported in 
previous studies fT71fT9l and is referred to as the 'hole-doping 
case' . It has been interpreted as resulting from pinning of the 
graphene hole-carrier density under the metal and the forma- 
tion of the surface-potential profile by charge transfer across 
the metal/graphene junction. Therefore, as V g is varied, the 
band profile changes from pp' p to pn'p, where p (n) and p f 
(n') refer to positive (negative) charge densities in the region 
under the metal segment or between the source and drain elec- 
trodes, respectively, as a function of the gate voltage. As a 
result, when V g < VDirac (Vg > Voirac), a pp'p (pn'p) junction 
is formed ifTTlfol 



Figure 1(d) shows how the conductance evolves in a mag- 
netic field as a function of V g , with V S( j = 100 pV and T = 1.7 
K. Near the Dirac point, conductance is strongly suppressed 
by the magnetic field. Furthermore, several resonance peaks 
emerge in the hole- and electron-carrier regions. In particu- 
lar, clear resonance peaks are observed in the hole-carrier re- 
gion (V g < Voirac)- These peaks correspond to the Coulomb- 
blockade effect, which is discussed below. The fluctuations 
in conductance as a function of V g are stable over time and 
reproducible. An enlarged view of the conductance as a func- 
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FIG. 2. (color). Coulomb blockade and quantum confinement under 
magnetic fields, (a) Large-scale current mapping as a function of V sd 
and V g at T = 0.23 K and B = 7.5 T. (b) High-resolution measure- 
ment of the differential conductance dI/dV sd as a function of V sd and 
V g for the green-boxed region of (a). Orange arrows point toward 
lines running parallel to the diamond edges, which correspond to the 
excitation spectra. The number of holes is fixed in each diamond 
region, and the '0, lh, 2h, • • • ' labels indicate the number of hole 
carriers in the first Landau level (also see the text), (c) Addition en- 
ergy E add as a function of hole number at B = 7.5 T. Inset shows the 
gate capacitance C g as a function of hole number at B = 7.5 T. (d) 
Magnetic field dependence of the addition energy, which is deduced 
from the size of the diamonds at each fixed magnetic field. 



tion of V g and B is shown in Fig. 3(a) for V S d = 500 pV and 
T = 0.23 K. 

To begin our discussion, we focus on the hole-carrier re- 
gion. Figure 2(a) shows a large-scale plot of the source- 
drain current (/) as a function of both V S d and V g at fixed 
B = 7.5 T and T = 0.23 K. Diamond- shaped regions appear 
in the hole-carrier regions, which indicates that the current is 
suppressed due to Coulomb blockade. Figure 2(b) shows a 
high-resolution measurement of the differential conductance 
dI/dV s & as a function of both V S d and V g that corresponds 
to the green-boxed region in Fig. 2(a). Characteristic fea- 
tures called 'Coulomb diamonds' manifest themselves, and 
the excitation spectra also appear, as shown by the orange 
arrows in Fig. 2(b). These data allow us to estimate the 
single-particle energy-level spacings (AE) to be « 2 meV. The 
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lines for the excited state are broader than thermal broaden- 
ing, which might be because of energy-dependent coupling of 
the excited state in the dot to the lead (7). 

Under a perpendicular magnetic field, a Landau level 
(LL) generally forms in the 2DEG system. The LL ener- 
gies (En) of the single-layer graphene sheet are given by 
E N = sgn(N)v? ^2eh\N\B, where e is the fundamental unit 
of charge, ft - h/(2n) where h is Planck's constant, is the 
Fermi velocity, and the integer N is the electron (N > 0) or 
a hole (N < 0) LL index EDJ. Note that the N = LL (E ) 
is independent of magnetic field in a single-layer graphene; 
i.e. Eq is fixed at the Dirac-point energy with the magnetic 
field. However, the N = LL spectrum does not appear in our 
device, which is attributed to the fluctuating spatial-disorder 
potential near the Dirac point and the pn'p band profile, as 
discussed below. In Fig. 2b, we show the Coulomb diamond 
associated with the ridge of the first LL of hole carriers, which 
is quantized into discrete single-particle levels. The '0, lh, 2h, 
• • • ' in each diamond indicates the number of carriers in the 
first LL of hole carriers. 

We now focus on the 'confinement-deconfinement transi- 
tion' controlled by the magnetic field. As discussed in Ref. 
I2TL such a transition between open-to-closed trajectories 
of Dirac particles is ubiquitous in this type of device struc- 
ture. Figs. 3(b) and (c) show the typical carrier trajecto- 
ries between contacts in graphene with a confinement poten- 
tial. For simplicity, we assumed that the electrostatic poten- 
tial between contacts formed by the metal/graphene junction 
is a parabolic function ED [221 : U = -(4U /L 2 )x 2 where 
Uq = Tivy ^CglVg - VDirad/^ (Cg is the gate capacitance per 
area), which corresponds to the potential depth at a fixed gate 
voltage. Here we consider that electron motion in the y di- 
rection is geometrically confined in our device. We set our 
model parameters on the basis of the measured area between 
the contacts of approximately 0.11 yum 2 (= 0.2 /zm x 0.55 
lim) and C % - 384 aF //mi 2 . The solution for trajectory was 
obtained by numerically integrating the classical equation for 
a Dirac particle in a magnetic field. The qualitative features 
of our picture are not sensitive to the details of confining po- 
tential within the window of Figs. 3(b) and (c), so far as the 
potential is smooth under a magnetic field (2TJI23). As shown 
in Fig. 3(c), all the trajectories are open for a sufficiently weak 
magnetic field, which corresponds to a continuum spectrum. 
As the magnetic field approaches the critical magnetic field 
(B c ), the orbits asymptotically approach the critical separatrix 
lines which correspond to the dashed black lines in Fig. 3(b). 
The orbits finally merge with the separatrix at B = B c . Above 
Z? c , closed orbits emerge that correspond to quasi-bound states 
and coexist with open trajectories, as shown in Fig. 3(b). At 
high magnetic fields, the closed orbits are approximately cir- 
cular. The quasi-bound states should be consistent with the 
Bohr-Sommerfeld condition (21] [22) and lead to a discrete 
energy spectrum. Furthermore, tunneling between open and 
closed trajectories leads to a finite resonance lifetime and QD 
formation (22l [23). Resonance shapes are governed by charg- 
ing effects via Coulomb blockade. Therefore, we consider 
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FIG. 3. (color). Field-induced quantum confinement-deconfinement 
in graphene. (a) Conductance (for V sd = 500 yt/V) mapping as a 
function of V g and B at T = 0.23 K. (b) and (c) Typical classical 
trajectories of Dirac particle in a parabolic confinement potential. At 
a high magnetic field (b), closed orbits are approximately circular and 
coexist with open trajectories. As the magnetic field approaches the 
critical magnetic field (B c ), the orbits become increasingly elongated 
until finally merging with the separatrix (dashed black lines) at B = 
B c . For a sufficiently weak magnetic field (c), all trajectories are 
open, which corresponds to a continuum spectrum. Based on the 
experimental data, we set the model parameters (see text) as follows: 
U = 0.2 eV (V g = -10 V), B = 1 (blue) and 2 T (red). 



that a genuine Dirac particle has been confined. In addition, 
the appearance of four clear Coulomb peaks, e.g. for the 7.5 
T [see the hole-carrier regions in Fig. 1(d)], suggests two- 
valley degenerate quantum states with two-fold spin degener- 
acy MM E3- 

Having understood the qualitative behavior of field-induced 
quantum dots in the graphene device, we quantitatively ana- 
lyze the Coulomb-blockade effect to substantiate our picture 
and better understand the finer details. The main panel of Fig. 
2(c) shows the addition energy (£ a dd) and gate capacitance of 
a single-electron charge (C g ) (inset) as a function of the num- 
ber of holes. The constant-interaction model of a QD gives 
E a dd = Eq+AE, where Eq = e 2 /C% is the charging energy and 
C% = C s +Cd+Cg, where C s and are the source and drain ca- 
pacitances, respectively [26]. Based on the experimental data 
shown in Fig. 2(b), we deduced £ a dd from the diamond size 
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and C g [= e/AV g (l + AE/E C )] from the gate- voltage spacing 
AV g between diamonds at V S( j = 0. By setting the energy-level 
spacings AE « 0, the classical capacitive-voltage relation for 
a single-electron charge becomes C g = e/AV g . From the exci- 
tation spectra of Fig. 2(b), AE is estimated for the lh and 2h 
Coulomb diamonds. Here £ a dd and C g for lh and 2h account 
for AE, and the other £ add and C g are deduced with AE « 0. 
Although the addition energy E a dd decreases rapidly as holes 
are added, C g is nearly constant or slightly increases. We find 
that C g (< 1 aF) decreases by two or three orders of magni- 
tude than the total capacitance (e.g. when Eq = e 2 /C^ = 4 
meV, the corresponding = 40 aF), which implies that the 
capacitance between the lead and the QD is the main contrib- 
utor to total capacitance. By using the C g data and the simple 
plate-capacitor model, the QD diameter is estimated to be 25- 
40 nm. These results are consistent with the above picture 
based on classical closed trajectory [Fig. 3(b)]. Figure 2(d) 
shows the magnetic field dependence of the addition energy, 
which is deduced from the size of diamond at each magnetic 
field value. The slope of each curve is ~ 1 meV/T; the slope 
cannot be attributed to the Zeeman effect, which would result 
in a slope of gpB = 0.12 meV/T, assuming a g factor of 2. 
The addition energy is mainly dominated by the charging en- 
ergy, which in turn is dominated by the capacitance between 
the lead and the QD. Therefore, these results indicate that cou- 
pling between the lead and the QD is sensitive to the magnetic 
field, which is consistent with our scenario discussed above. 

Finally, we consider the electron-carrier regions (V g > 
^Dirac), in which the conductance peaks are irregular functions 
of the magnetic field [Fig. 1(d)]. In particular, peak crossings 
occur in Fig. 3(a), which implies transport through a multi- 
QD system. We consider that a multi-QD system forms as a 
result of the inhomogeneous pn'p band profile along the con- 
tact metals (in the y direction) and the spatial disorder poten- 
tial fluctuations that arise because of the electron-hole puddle 
near the Dirac point (27l |28) . The formation of such a multi- 
QD system can also be deduced from the transport character- 
istics at a fixed magnetic field, as shown in Fig. 2(a). We 
observe a large transport gap 20 meV) between the Dirac 
point and the lh states in the ridge of the first LL, but the 
diamond structure collapses near the Dirac-point side and an 
irregular saw-tooth pattern is observed in the electron-carrier 
region. 

In conclusion, we propose a graphene device structure that 
is not disturbed by the structural fluctuations of constric- 
tions and that is directly connected to the mesoscopic 2DEG 
system. We demonstrated the Coulomb-blockade effect and 
quantum confinement, which are induced by both a uniform 
magnetic field perpendicular to the graphene sheet and an 
electrostatic surface-potential formed by the metal/graphene 
junction. Our experimental results indicate that a Dirac- 
type particle is confined and that quantum confinement- 
deconfinement transition occurs in the mesoscopic graphene 
system. These findings open the door to the realization of 



a single-Dirac-fermion device and the demonstration of rela- 
tivistic quantum physics in a condensed-matter system. 
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